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WHAT IS CLAIMED IS: 



1. A charge rate estimating apparatus for a 

secondary cell, comprising: 

a current detecting section capable of measuring 
a current flowing through the secondary cell; 

a terminal voltage detecting section capable of 
measuring a voltage across terminals of the secondary 
cell ; 

a parameter estimating section that calculates 
an adaptive digital filtering using a cell model in a 
continuous time series shown in an equation (1) 
estimates all of parameters at one time, the 
parameters corresponding to an open-circuit voltage 
which is an offset term of the equation (1) and 
coefficients of A(s), B(s). and C(s) which are 
transient terms; and 

a charge rate estimating section that estimates 
the charge rate from a relationship between a 
previously derived open-circuit voltage V 0 and the 
charge rate SOC using the open-circuit voltage V 0 . 

v = »V (1). wherein s denotes a 

A(s) C(s) ° 

Laplace transform operator, A(s), B(s), and C(s) 
denote poly-nominal functions of s. 

2. A charge rate estimating apparatus for a 

secondary cell as claimed in claim 1, wherein the 
open-circuit voltage V 0 of the cell model in the 
continuous time series shown in the equation (1) is 
approximated by means of an equation (2) to provide 
an equation (3) and the digital filter calculation is 
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carried out using the equation ( 3 ) and equivalent 
equation (4), h is estimated in at least equation (4), 
the estimated h is substituted into equation (2) to 
derive an open-circuit voltage V 0 , and the charge 
rate is estimated from a relationship between the 
previously derived open-circuit voltage V 0 , and the 
charge rate is estimated from a relationship between 
the previously proposed open-circuit voltage V 0 and 
the charge rate (SOC). 



Vo = --- (2) 

s 



v = ( B & , 1 -*W = s*B{s)'C(s) + h*A(s) mI _ 
1 A(s) C(s) sj s • A(s) • C(s) 



s*A(s)'C(s) mV = s*B(s)*C(s)+h'A(s) mI _ 

G x (s) G x (s) 



wherein s denotes the Laplace transform operator, 
A(s), B(s), and C(s) denote poly-nominal equation 
functions, h denotes a variable, and 1/Gi(s) denotes 
a transfer function having a low pass filter 
20 characteristic. 



3. A charge rate estimating apparatus for a 
secondary cell as claimed in claim 1, wherein the 
open-circuit voltage V 0 of the cell model in the time 

25 continuous time series is approximated in an equation 
(2) to calculate an equation (3), the adaptive 
digital filter calculation is carried out using an 
equation (4) which is equivalent to the equation (3), 
A(s), B(s), and C(s) are estimated from equation (4), 

30 the estimated A(s), B(s), and C(s) are substituted 



-25- 



into equation (5) to determine V 0 /G 2 (s) and the 
charge rate is estimated from the relationship 
between the previously derived open-circuit voltage 
V 0 and the charge rate (SOC) using the derived 
5 V 0 /G 2 (s) in place of the open-circuit voltage V 0 , 



V 0 = -•/ --- (2) 

s 



10 



V - 



(B(s) f 1 ^\ _ S *B(s)'C(s) + h'A(s) 
A(s) C(5) s) s*A(s)*C(s) 



s-A(s)-C(s) _ v = s*B(s)*C(s}+h'A(s) 



Vr 



o_ = C ( 5 ) .Qf_i(?l.j) (5)/ wherein s denotes the 



G 2 (s) G 2 (s) A(s) 
Laplace transform operator, A(s) , B(s), and C(s) 
15 denote the poly-nominal (equation) function of s, h 
denotes a variable, 1/Gi(s) and 1/G 2 (s) denote 
transfer functions having the low pass filter 
char a c t eris t ic s 

20 4 . A charge rate estimating apparatus for a 

secondary cell as claimed in claim 1, wherein the 
cell model is calculated from an equation (6), 

K'(T, •s+l) T 1 

V = = •/ + V 0 , wherein K denotes an 

T x • s + 1 T 3 • s + 1 

internal resistance of the secondary cell, T i# T 2 , 
25 and T 3 denote time constants, 1/Gi(s) denotes a low 

pass filter having a third order or more, and 1/G 2 (s) 
denotes another low pass filter having a second order 
or more . 
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5. A charge rate estimating apparatus for a 
secondary cell as claimed in claim 4, wherein A(s) 
T x • s + 1 # B(s) = K • (T 2 + 1), C(s) - T 3 • s + 1. 

6. A charge rate estimating apparatus for a 
secondary cell as claimed in claim 5, wherein 

V - *•<?•. " + '> . / . ±„ ....,„ 



(a • s 3 + b • s 2 +s) • V = (c • s 3 + d • s 2 
+ f ) •! ( 10) . 



7. A charge rate estimating apparatus for a 
secondary cell as claimed in claim 6 , wherein a = T 
• T 3 . b = Ti + T 3 , c = K • T 2 • T 3 , d = K • ( T 2 + T 3 ), e 
= K + A • T 1# f = A (11). 

8. A charge rate estimating apparatus for a 
secondary cell as claimed in claim 7, wherein a 
stable low pass filter G x (s) is introduced into both 
sides of the equation (10) to derive the following 
equation : 

— ^— (a-s 3 + b*s 2 + s)*V = — 4t(c*s 3 + d*s 2 + e*s + f)*i 
(12) . 



9. A charge rate estimating apparatus for a 
Secondary cell as claimed in claim 8 , wherein 
actually measurable currents I and terminal voltages 
V which are processed by means of a low pass filter 
are as follows : 
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Ii = 



I 2 = 



I3 = 



G,(s) 

s 2 
G,(s) 

3 



G,(s) 

1 



Vi = 



V, = 



v 3 = 



(7,(5) 



V , 

'V . and 



- (13). 



15 



10. A charge rate estimating apparatus for a 
secondary cell as claimed in claim 9, wherein, using 
the equation (13), the equation of (12) is rewritten 
and rearranged as follows : 



10 Vi = [V; 



Iol = 



(15) and 



- a 
-b 

c 
d 
e 

f 

the equation (15) corresponds to a general equation 
which is coincident with a standard form of a general 
adaptive digital filter of equation (16): y = co T • 6 

(16), wherein y = Vi , co T = [V 3 V 2 I 3 I 2 Ii I 0 J. 

and 

-a 
-b 

c 
d 
e 

f 



e = 



--- ( 17) . 



ii . 



A charge rate estimating apparatus for a 
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Secondary cell as claimed in claim 10, wherein a 
parameter estimating algorithm with the equation (16) 
as a prerequisite is defined as follows: 



Y(k) 
6(k) 
P(k) = 

Xi(k) = 



A 3 (*) 



1+ A 3 (k)»co T (k) • P(k - 1) • co(k) 

6(k -1) - Y(k)-P(k-l)*w(k)»[co T (k)« 6(k -l-y(k)] 



P(yt-1)-^ 



A 3 (k) • - 1) • co(k) •a) T (k)» P(k - 1) 



trace\P'(k) } 



1 + A 3 (A:) F(A: - 1) • co(k) 

trace\p'(k) } 



A, (*) 



Yu 



trace 



trace 



jf (*) } 



{ 



Yu 



trace 



\P\fc) } 



trace 



{*»'(*) } 



10 (18), wherein 0(k) denotes a parameter estimated 

value at a time point of k (k = 0, 1, 2, 3 ), Xi , 

X 3 (k), yu , and y l denote initial set value, h < X 1 < 1 , 
0 < X. 3 (k) < oo , P(0) is a sufficiently large value, 
0(0) provides an initial value which is non-zero but 

15 very sufficiently small value, and tracejp} means a 
trace of matrix P. 



12. A charge rate estimating method for a secondary 
cell, comprising: 
20 measuring a current flowing through the 

secondary cell; 

measuring a voltage across terminals of the 
secondary cell; 

calculating an adaptive digital filtering using 
25 a cell model in a continuous time series shown in an 
equation ( 1 ) ; 



10 
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estimating all of parameters at one time, the 
parameters corresponding to an open-circuit voltage 
which is an offset term of the equation (1) and 
coefficients of A(s), B(s) # and C(s) which are 
transient terms; and 

estimating the charge rate from a relationship 
between a previously derived open-circuit voltage V 0 
and the charge rate SOC using the open-circuit 
voltage V 0 , 



V = • /+ 1 m V n (1), wherein s denotes a 

A(s) C(5) 

Laplace transform operator, A(s), B(s), and C(s) 
denote poly-nominal functions of s. 



15 13. A charge rate estimating method for a secondary 
cell as claimed in claim 12, wherein the open- 
circuit voltage V 0 of the cell model in the 
continuous time series shown in the equation (1) is 
approximated by means of an equation (2) to provide 

20 an equation (3) and the digital filter calculation is 
carried out using the equation (3) and equivalent 
equation (4), h is estimated in at least equation (4), 
the estimated value of h is substituted into equation 
(2) to derive an open-circuit voltage V 0 , and the 

25 charge rate is estimated from a relationship between 
the previously derived open-circuit voltage V 0 , and 
the charge rate is estimated from a relationship 
between the previously proposed open-circuit voltage 
V 0 and the charge rate (SOC). 

30 



(2) 
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V = 



'*£>+_!_. *W - s-B(s)-C( S ) + h-A(s) 

A(S) C(5) 5 J 5^(5)-C(5) 



s'A(s)-C(s) _ v = s-2?(*)-C(s>f-/»-,4(s) 

wherein s denotes the Laplace transform operator, 
A(s) # B(s), and C(s) denote poly-nominal equation 
functions, h denotes a variable, and 1/Gi(s) denotes 
a transfer function having a low pass filter 
characteristic . 



10 14. A charge rate estimating method for a secondary 
cell as claimed in claim 12 , wherein the open- 
circuit voltage V 0 of the cell model in the 
continuous time series is approximated in an equation 
(2) to calculate an equation (3), the adaptive 

15 digital filter calculation is carried out using an 

equation (4) which is equivalent to the equation (3), 
A(s), B(s), and C(s) are estimated from the equation 
(4), the estimated A(s), B(s), and C(s) are 
substituted into equation (5) to determine V 0 /G 2 (s) 

20 and the charge rate is estimated from the 

relationship between the previously derived open- 
circuit voltage V 0 and the charge rate (SOC) using 
the derived V 0 /G 2 (s) in place of the open-circuit 
voltage V 0 , 

25 

Vo = -•/ --- (2) 



V = 



f B(s) 
Ms) C(. 



(5) S) 



s • A(s) • C(s) 



(3) 
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5 * A(s) • C(S) 

G t (s) 



•V 



s*B(s)*C(s)+h'A(s) 



(4) . 




wherein s denotes the 



Laplace transform operator, A(s), B(s), and C(s) 
denote the poly- nominal (equation) function of s, h 
denotes a variable, 1/Gi(s) and 1/G 2 (s) denote 
transfer functions having the low pass filter 
characteristics 

15. A charge rate estimating method for a 
secondary cell as claimed in claim 12, wherein the 
cell model is calculated from an equation (6), 

K*(T, • j + r 1 
V = — '.•/ + V 0 , wherein K denotes an 

^•5+1 r 3 **+i 
internal resistance of the secondary cell, T i# T 2 , 
and T 3 denote time constants, 1/Gi(s) denotes a low 
pass filter having a third order or more, and 1/G 2 (s) 
denotes another low pass filter having a second order 
or more . 

16, A charge rate estimating method for a secondary 
cell, comprising: 

measuring a current I(k) flowing through the 
secondary cell; 

measuring a terminal voltage V(k) across the 
secondary cell; 

storing the terminal voltage V(k) when a current 
is zeroed as an initial value of the terminal voltage 
AV(k) = V(k) - V_ini; 
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20 



determining instantaneous current values I 0 (k), 
Ii(k), and I 3 (k) and instantaneous terminal voltages 
Vi(k), V 2 (k), and V 3 ( k ) from an equation (19), 

1 



In = 



II = 



I 2 = 



I 3 = 



S 

s~ 



'I . 

I . 



0,(5) 

s 3 

0,(5) 



Vi = 



V 2 = 



v 3 = 



0,(5) 



V , 



G,(5) 



•V 



G,(s) 



'V , and 



- (19), wherein pi denotes a 



o,(s) (pi-s+iy 

constant determining a responsive characteristic of 
10 G!(s); 

substituting the instantaneous current values 
I 0 (k), I x (k), and I 3 (k) and the instantaneous terminal 
voltages Vi(k)» V 2 ( k ) , and V 3 (k) into an equation (18), 



Y(k) 
15 8(k) 



\ 3 (k) 



1 + A 3 (k) • co 7 (k) • P(k - 1) • co(k) 



= 6(k -1) - Y(k)*P(k-l)'to(k)»[o3 T (k)«0(k -l) -y(k)] 



P(k) = 



Xi(k) = 



\(k) 
trace 



p(k-i) 



A 3 (k) ' P(k - 1) « co(k) • co T (k) « P(k - 1) 
1 + A 3 (k) • co 7 (k) * P(k - 1) • co(k) 



P\k) 
A,(k) 



km } 



Yu 



trace 



\P'(k) } 



Yu 



{ ^1 



trace 



\P\k) } 



trace 



fr'(*) I 



Yu Yi 

| trace\p\k) } . trace\p'(k) } ^ ^ 
Yl Yl 

(18), wherein 6(k) denotes a parameter estimated 

value at a time point of k (k = 0, 1, 2, 3 ), Xi , 

^3 ( k ) , y u ' and yL denote initial set value, b < Xi < 1 , 
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0 < X 3 (k) < oo. P(0) is a sufficiently large value, 
0(0) provides an initial value which is non-zero but 
very sufficiently small value, trace{P} means a trace 
of matrix P, wherein y(k) = V x (k) 



co T (k) = [V 3 (k) V 2 (k) I 3 (k) I 2 (k) Ix(k) 



Io(k) ] 



8(k) = 



(20); 



-a{k)] 
-b(k) 
c(k) 
d(k) 
e(k) 
/(*) 

substituting a, b, c, d , e, and f in the 
parameter estimated value 6(k) into and equation (22) 
to calculate V 0 ' which is an alternate of V 0 which 
corresponds to a variation AV 0 ( k ) of the open-circuit 
voltage estimated value from a time at which the 
estimated calculation start is carried out; 



Vo 



+ b 



+ V 4 -c 



- d 



G 2 (s) 

e • I 4 ( 22 ) ; and 

calculating an open-circuit voltage estimated 
value V 0 (k) according the variation AV 0 ( k ) of the 
open-circuit voltage estimated value and the 
terminal voltage initial value V_ini, 

17. A charge rate estimating method for a secondary 
cell as claimed in claim 16, which further comprises: 

calculating a charge rate (SOC) from the open- 
circuit voltage estimated value AV 0 ( k ) . 



25 
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18. A charge rate estimating method for a secondary 
cell as claimed in claim 17, wherein the charge rate 
(SOC) is calculated using a correlation map between 
the open-circuit voltage V 0 and the charge rate of 
the secondary cell from the open-circuit voltage 
estimated value AV 0 ( k ) . 
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